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We point out that two different definitions of the super- 
fluid density - through statistical response to static gauge 
phase and through dynamic response to altering gauge phase 
- yield, generally speaking, different quantities in d < 3. The 
physics leading to this difference is associated with the equi- 
librium statistics of supercurrent states. Some experimentally 
observable consequences of this fact are discussed. 

PACS numbers: 67.40.-w, 67.40.Rp, 74.76.-w 

In this report we would like to discuss the relation be- 
tween two alternative definitions of the superfluid den- 
sity. Consider a system with ring geometry. One way 
to define the superfluid density is to consider system re- 
sponse to the gauge field generated by the infinitesimal 
flux 4> threading the ring (we use units such that h = 1 
and particle mass m = 1) 
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where T is the free-energy, L is the system size, d is the 
system dimensionality (for simplicity we assume that the 
system size is the same in all directions - other examples 
are discussed at the end of the paper), and ip = 2tt(4>/ <f>a) 
, where (j)Q is the flux quantum. We do not specify here 
what is the value of the flux quantum, since the gauge 
field considered is not necessarily of elecrtro- magnetic ori- 
gin. For neutral systems (helium, ultracold atomic gases) 
it can be introduced by rotating the system. More gener- 
ally, we may ascribe to the particles fictitious charge and 
couple it to fictitious gauge field - the idea is to look at 
the response of the system to the twisted boundary condi- 
tions on the ring (we choose axis x as the direction around 



the ring): \E'(. 
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[here Xj is the ^-coordinate of the j-th particle]. In fact, 
this relation may be considered as our initial definition 
of the phase tp. The equilibrium particle number current 
is given then by 



= p<P v /L 
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On another hand, one may define superfluid density 
as a coefficient in effective.-long-wavelength action which 
governs phase fluctuations^: 



Ps_ 
2 



dV (V$) 2 



(3) 



We do not intend to discuss here under what conditions it 
is possible to introduce phase field in the long-wavelength 
limit - this is rather well known and studied problem (see, 



for example, Refs. |l],|2|). In most general terms, Eq. (|3j) 
works when V$ is a well defined operator. In dimensions 
d > 1 this would require the existence of the topolog- 
ical long-range order. In ID, where the concept of the 
long-range order at any finite temperature is meaning- 
less, Eq. (|^) is valid provided the low-energy spectrum 
of the system is gapless (metallic state) and relaxation 
times of current states, 1/r, are very long. In particular, 
if we label current states in ID by index /, where 
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then the sufficient condition reads 1/tj <C 2irc/ (LK), 
where c is the sound velocity, and K/2 is the one-particle 
density matrix index, p(x) ~ x~ K l 2 at zero temperature. 
As we discuss below this condition is satisfied at low tem- 
peratures, including temperatures exceeding finite-size 
quantization limit T 3> 2nc/L. 

Both definitions are used in literature with the unques- 
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tioned assumption that Ps Ps- ^ n wnaTi follows, we 
demonstrate that although in 3D one indeed may, with 

(F) 

macroscopic accuracy, make no distinction between p s 
and p s , in 2D and especially in ID, at finite tempera- 

(F) 

ture p s 7^ Ps- We then discuss some interesting con- 
sequences associated with this inequality, and how these 
may be directly tested experimentally. 
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We start from the well-known relation between p s 
and the statistics of worldline winding numbers, M, (see 
Ref. |) 
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where M is integer. Formally, Eq. (g) immediately fol- 
lows from the 27r-periodicity in ip of the partition func- 
tion Z(jp), implying that it can be Fourier-expanded as 
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The statistical average in 



Eq. (ph is taken with the distribution function Zm- The 
actual extreme usefulness of these formulae for numeri- 
cal simulations comes from the fact that the numbers M 
are directly related to the topology of the closed particle 
worldlines in imaginary timefl 

On the other hand, the partition function Z(ip) can be 
calculated directly from the effective long-range action 
([|), employing gauge- invariance argument that external 
vector-potential and/or twisted boundary condition may 
be transformed into the phase of the wavefunction. We 
thus write 



Z(ip)(x Jv$ exp{-F cff [$]} 
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where $ is parameterized as 
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We then immediately see that 
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where Zq is independent of ip. 

The second derivative of —TlnZ(ip) at <p = 0, related 
to PgL d ~ 2 , then yields 
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the averaging being performed with the normalized dis- 
tribution 
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Eqs. ( |TT| ) and (|T^) establish a relation between p^ and 
p s . One can also directly relate p s to the statistics of 
worldlinc winding numbers. Using the known properties 
of 03-functions, the following identity is valid 
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where M's are some integers, which, by the uniqueness 
of the Fourier transform, are immediately identified with 
the worldline winding numbers. We thus see that the 
worldline winding numbers have Gaussian distribution 
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As is clear from Eq. (|ll| ), pg and p s in general dif- 
fer. However, at T — > and finite 7, i.e. in the "meso- 

(F} 

scopic" limit, we find that p s — p s . Moreover, in 3D 
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case, one may safely replace p s with p s since, accord- 
ing to Eqs. (jll]-(lf), finite -size corrections are exponen- 
tially small in macroscopic parameter: the probability of 
finding non-zero I = ±1 dies away as exp{— 2tt 2 p s L/T}. 
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This is probably the reason why it became a custom to 

{ F) 

identify p s and p s . 

In ID systems the situation changes drastically. In 
the thermodynamic limit, when L — > oc first, that is 
TL — > oo at any finite temperature, we have negligible 
p { s F) [~ LTexp{- J LT/2p s }, see Eq. Q)] while p s is fi- 
nite. At this point we would like to emphasize that zero 

{ F) 

p s is not in contradiction with the fact that the system 
may remain superfluid in any dynamic sense. Equations 
( [TT[]T^ ) make it clear that p^ — > simply due to the 
broad distribution over the current states, with typical 
I ~ (LT) 1 / 2 ^> 1. However, |i,elaxation times of current 
states in ID may be very longp (see also below), up to as- 
tronomical scale. Therefore, out-of-equilibrium state ob- 
tained experimentally, say, by switching on/off the gauge 
phase ip, will support a supercurrent j ~ p s ip/L. 

To estimate relaxation times tj in ID rings (obviously, 
for the relaxation time to be finite the translational in- 
variance must be violated by impurities, or lattice po- 
tential) one may use an effective Hamiltonian describ- 
ing low energy dynamics of supercurrent states and den- 
sity fluctuationso. Relaxation rate is obtained then from 
the golden-rule expression for the terms breaking trans- 
lational invariance. Formally, this problem is identical 
to the study of particle dynamics rwith Ohmic coupling 
to the oscillator bath environmental - an analogy widely 
used in ID systems since Ref. [| The standard result for 
the transition rate from 7 to / + 1 due to single impurity 
backscattering potential, A, is given by (see, e.g., Ref. u\) 
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where oj is some model-dependent high-energy cut-off. 
If parameter K is small, then at low temperature one 
may easily get extremely long tj. The answer for the 
disordered ring with white-noise disorder is just the single 
impurity expression ( |l5| ) multiplied by L. 

One may directly test the striking difference between 

( F) 

p s and p s in ID systems experimentally. In fact, we 
propose a unique measuring device for magnetic fields, 
free of limitations set by quantization of magnetic flux. 
Suppose we would like to measure local magnetic field at 
the point R by placing at this point a device recording 
magnetic flux through some small surface with subse- 
quent reading the information recorded. We point out 
that such a device can be realized as an ensemble of ID 
(quasi-)superconducting rings. Being placed at the point 
R at high temperature (say, T ~ u>o) when relaxation 
times are short, and then cooled down to low temper- 
atures when relaxation times become very long (much 
longer than expected "storage time"), the ensemble will 
not respond to the local field, since in this case the re- 

( F) 

sponse is statistical and p s ~ 0. However, removing 
cooled device from the point R we excite finite super- 
current, since now we probe dynamic response p s . It is 
important to note, that, in contrast to the case of mas- 
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sive superconducting ring, the stored magnetic flux is not 
quantized. 

Another interesting experiment directly follows from 
Eq. (|l^). A single ID ring being cooled down to low 
temperature will typically end up in a state with non-zero 
supercurrent j ~ (p s T c ^/ L) 1 / 2 , where T cx is the "freez- 
ing" temperature at which the current relaxation time 
becomes comparable to the experimental cooling time. 

The 2D case is even more intriguing. For an isotropic 
system with L x = L y , Eqs. (fLl)), (|||) and (jT|) ar e 
system-size independent, in particular 
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The famous Nelson-Kosterlitz formula 

7T 



( F) 

relates transition temperature and p s [not p s ! 
Ref. § and also Ref.§). Substituting it into Eq. (|l 
find that with very high accuracy 



p { s F, (T c )=p s (l~16ne^) 
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(see 
), we 
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For purely numerical reasons the relative difference 
{p^ — Ps)l Ps < ^ • 10~ 4 is very small. That is why in a 
recent state-of-the-art numerical study of the Berezinskii- 
Kosterlitz-Thouless transitionEI formally incorrect rela- 
tion T c = (tt/2)p ( ;P (or (M 2 ) = 2/ttB) was successfully 
used to fit the data. However, unlike in 3D case, the dif- 

ference between p s and p s is not controlled any more 
by small parameters. Furthermore, if a 2D system is 
elongated geometrically: L x 



> L y , then 
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where index x for p s and Wj is used to specify that we 
are dealing with the response to the twisted boundary 
condition along ir-axis. Clearly, at L x 3> L y the physics 
becomes analogous to that of ID rings. We thus see that 
in 2D pg depends on the sample geometry and at L x > 

is not related to T c , Eq. (17), 



4tt£„ 



the value of pg 



even approximately. 

It seems rather easy experimentally to confirm these 
conclusions by producing narrow Carbino-disk samples. 
Since in this case we are dealing with macroscopic phase 
transition phenomenon, at T < T c the current relax- 
ation times become astronomically long and the above- 
mentioned experiments on ID rings would simply require 
cooling below T c . 
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